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Resonant  Frequencies  and  Damping  of  a  Liquid  Drop  With  a  Rigid  Spherical  Core 


by  C.  A.  Morrison 
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Abstract 


Rayleigh's  method  is  used  to  obtain  the  resonant  frequencies  of  a  liquid 
drop  with  a  rigid  core.  The  decay  time  of  the  drop  is  calculated  by  using 
Rayleigh's  dissipation  function  and  by  assuming  that  the  fluid  in  the  drop  and 
the  external  medium  are  viscous.  It  is  shown  that  special  cases  reduce  to 
known  results. 
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In  the  formation  of  small  hall  within  the  atmosphere,  conditions  exist  in 
which  the  inner  core  is  ice,  with  a  layer  of  water  on  the  outside.  Such  a 
composite  drop  has  resonant  frequencies,  and  if  both  the  air  surrounding  the 
drop  and  the  liquid  of  the  drop  are  viscous,  these  resonances  are  damped. 

Such  resonances  and  perhaps  the  damping  may  be  observed  by  techniques  used  to 
observe  oscillations  in  completely  liquid  drops.1**3 

In  tills  report  we  investigate  the  resonant  frequencies  and  damping  for  a 
composite  drop.  Both  the  liquid  of  the  drop  and  the  external  medium  (air) 
will  be  assumed  to  be  viscous.  The  technique  used  is  *hat  of  Rayleigh,4  which 
we  have  used  previously  for  determining  the  resonant  frequencies  of  a  liquid 
drop  in  an  electric  field.5  The  damping  will  be  treated  in  a  manner  similar 
to  that  used  in  Morrison  et  al.5 

2.  KINETIC  AND  POTENTIAL  ENERGY  OP  THE  DROP 

The  inner  rigid  core  of  radius  b  is  surrounded  by  a  fluid  of 

density  p^  and  viscosity  The  exterior  of  the  drop  of  quiescent  radius  a 

is  surrounded  by  fluid  or  air  of  density  p  and  viscosity  n  .  As  in  our 

o  o 

previous  analysis,5  we  assume  that  the  outer  surface  of  the  drop  is  givenby* 

r(0,t)  »  a  (t)  ♦  T  a  (t)  P  (cos  a)  ,  (2.11) 

o  n  n 

n 


^Equation  numbers  containing  an  I  are  found  by  that  number  in  reference  (5). 
Thus,  12.11)  is  (2.1)  in  reference  (5). 
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where  the  prime  on  the  sum  indicates  that  the  term  n  •  0  is  missing,  the 
constraint  that  the  volume  of  the  liquid  contained  in  the  outer  shell  is 
constant  is  the  same  as  in  Morrison  et  al. 5  To  show  this,  we  see  that  the 
volume  of  the  liquid  in  the  drop  is 


1  H3 

v  -  2*  /  -  % )  dM  , 

-1  3  3 


where  d|i  -  -d(cos  0),  and  the  equilibrium  volume  is  given  by 


V  <■  —j  (a3  -  b3)  . 


If  the  result  in  equation  (1)  is  equated  to  that  of  equation  (2),  the  volume 


of  the  inner  sphere  cancels,  and  the  constraint 


identical  to  that  of  a 


liquid  sphere,  giving 


•o  *  •C1 


-  —  l'  - )  . 

.2  L  2n  ♦  V  * 


(2.21) 


a  n 


The  energy  due  to  surfaee  interface  tension,  y,  at  the  outer  surface  is 
obviously  the  same  as  in  Morrison  et  al,s  or 


,2  *•  (a  -  IHn  ♦  2)  2* 

I,  •  «T  [.  +  l  2(to  ♦  *.)  • 


(2.31) 
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We  shall  assume  that  the  velocity  of  the  field  In  both  regions  Is 
derivable  free  a  velocity  potential#  so  that 


Tt  -  -  v*i 


▼  •  -  V# 

o  o 


b  <  r  <  r(0,t) 


r(0,t)  <  r  , 


(3) 


with  and  0Q  given  by  suitable  solutions  to  Laplace's  equation*  The  total 
kinetic  energy  Is  then  given  by 


1  -  7  Pt  /  (»♦*)*  4t  ♦  i  p„  /  (»*0)2dT  , 


(4) 


where  the  volume  Integrals  cover  the  regions  given  In  equation  (3).  The 
voluee  integrals  In  equation  (4)  can  be  converted  to  surface  Integrals  by 

/  (V#)2dT  -  /  ♦(?♦)  •  do  ,  (5) 


idiere  do  is  the  incremental  surface  area  directed  outward  from  the  volume* 
Oeing  this  result  in  the  first  Integral  in  equation  (4),  we  get 

/  at  -  2. 


/  dv  r  (0,t)  g. 


-1 


(4) 


3 


where  all  the  functions  in  the  integrand  are  evaluated  at  r(  0,  t)  and  where  we 
have  uaed  the  condition  »  0  at  t  ■  b.  The  second  integral  in  equation  (4) 
can  be  evaluated  in  the  sane  Banner  to  give 


3  1  3 

/  (V*  P  dt  -  -2*  /  dp  r  <B,t>  0 
-1 


o  0r 


(0»t) 


(7) 


where  we  have  assueed  that  *  0  at  r  ♦  ». 


The  results  given  in  aquations  (6)  and  (7)  can  be  used  to  obtain  the 
kinetic  energy  in  equation  (4)  once  the  velocity  potentials  are  determined. 
The  appropriate  solutions  of  Laplace's  equation  for  the  two  regions  are 


A 

♦i  ■  l  (“Jjjf  ♦  rnBn)  pn(co*  ,  b  <  r  <  r(0,t)  , 

n  t 

c 

♦  •  l  P  (cos  •)  •  r(0,t>  <  r  <  •  , 

n  rn 

(8) 

where  An,  Br,  and  CR  are  constants  to  be  determined  by  the  boundary 
conditions.  Also,  in  selecting  the  appropriate  solution  for  the  exterior 
region,  (r  >  r(  0,0),  we  have  chosen  the  solution  that  vanishes  at  infinity 
because  ms  assume  that  the  fluid  is  at  rest  far  from  the  drop. 
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Before  obtaining  explicit  expressions  for  the  constants  in  equation  (8), 
it  is  convenient  to  evaluate  the  tens  in  the  kinetic  energy  by  using  these 
expressions.  For  the  exterior  region,  we  use  the  second  equation  (8)  in 
equation  (7)  to  obtain 


J  (^0)2  dt 


-4*  l 

n 


In  »  1)  Cn 
(2n  +  1 )  2n+1  ' 

A 


(9) 


where  we  haws  taken  r  ■  a  in  the  integrand  of  equation  (7)»  higher  order  terns 

give  corrections  of  the  order  a3,  which  we  are  ignoring.  Using  equation  (8) 

n 

in  equation  (6),  w»  get  for  the  interior  region 


/  IT 


2n+1 

2n  +  1  t^n 


(n  ♦  1 )A  A 

2n+1  1  ^n  *  2n  +  0 

a  a 


(10) 


where,  as  in  equation  (9)  we  have  taken  r(0,t)  *  a  for  the  sane  reason  as 
given  in  the  derivation  of  equation  (9).  the  results  given  in  equations  (9) 
and  (10)  can  be  used  to  calculate  the  kinetic  energy  once  the  constants  are 
evaluated. 

Since  we  are  assuming  the  inner  sphere  to  be  rigid,  we  have  vf  ■  0  at 


r  ■  b,  and  by  using  equations  (3)  and  (8)  we  obtain 


w.'-L'-lrcs 
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/  (*♦„)  dt 


■  4*a3  l 


(n  +  1 )(2n  ♦  1 ) 


Finally,  the  kinetic  energy  is  given  by  using  equations  (14)  and  (15)  in 
equation  (4),  or 


i 


»  -  a..3  T  ~  “V  + _ Sa _ 

^  ‘•Mn  +1)(2n+1)F  (n  +  1)(2n  +  D^t 


In  tiie  absence  of  losses,  the  results  given  in  equation  (16)  and  in 


i 


equation  (2.31)  can  be  used  to  obtain  the  equation  of  notion  for  the  aR[l>  »  T 

-  0  and  “r  *  o]i  •••  equation  (2.81). 

dt  a;  3*n 
n 

•• 

_  p. (2n  +  1  -  nr  )  a 

•  P - Sf - ~  +  po)  Tr^TT  *  T<“  -’><”*  s,»„  -  0  •  <”> 


The  resonant  frequency,  ain,  for  aR  is  then  given  by 


y(n  -  1)(n  l)(n  +  2)nT 


n  3 


a  [pi(2n  ♦  1  -  nrn)  ♦  PonrJ 


Several  checks  on  the  validity  of  equation  (18)  can  now  be  nade.  If  in 

equation  (18)  we  let  b  *  0 ( r  ■  1)  and  p  -  0,  we  obtain  Rayleigh's  result,' 

n  o 


2  v(n  -  1 )n(n  »  2) 


V  %  ‘T’.  ••  .% 
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V  -*  -‘ 


1 


HDL-TL-83-3,  February  1 983 


Of  course,  if  b  -  a(T  >0)  and  there  is  no  fluid  to  oscillate,  consequently 

n 

u  -  0  results.  If  we  let  b  ■  0  in  equation  (18),  then  we  obtain  the  resonant 
n 

frequency  of  a  bubble  in  a  liquid  or  air  as 


-  1 )n(n  +  1 ) 


natn 


a  [<n  +  1)pA  ♦  npQ] 


and  this  is  identical  to  the  result  given  by  Laab.6 


3.  LOSSES  DOE  10  VISCOSITY 

As  in  Morrison  et  al5  we  shall  calculate  the  Rayleigh  dissipation  factor, 
R,  to  include  losses  due  to  viscosity  in  the  equation  of  notion  for  the 
an(t).  The  fore  given  by  Landau  and  Lif shits7  for  the  dissipation  factor  is 


R  n  /  (Vv2)  •  do 


where  n  is  the  viscosity  of  a  particular  region  and  do  is  the  increaental  area 
directed  outward  fron  the  volune  enclosed  by  the  surface. 

We  shall  consider  the  region  b  <  r  <  r(0,t)  first.  For  the  surface 
r  *  b,  we  have 


;-v-' 


%.ii  ^sJ*  -v  -.  V-J*  sAJ 


s5l 


HDL-TL-83-3,  February  1983 


m 


da  -  -2wrr  dy  , 


and  since  »  ■  0  at  r  «  b,  we  have 


2  2  8V0 
da  •  Vv  ■  -2wr  dy  —jj£  • 


Using  the  result  of  equation  (23)  in  equation  (21),  we  have 


,  (2n  +  1)a2 

/  (Tv2)  •  da  -  8w  l  - %_-r"r 

n  (n  ♦  1  Jna2"  2r2 
n 


for  the  surface  r  -  b.  For  the  surface  r  »  r(0,t),  we  have 


k 


/  (^) 


e*  I  — |  [» 

n  nr 
n 


-  ,  -  r1  . 

n  -f  1  »• 


w-..  ’’fr  1 >  r*>.|; 
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The  total  dissipation  factor 
combining  equations  (24)  and 


for  the  region  b  <  r  <  r(8,t)  is  then  given  by 
(25)  to  give 


•2 

4»V  £  ^  tn 

1  n  nrz 
n 


_  n(n  ♦  2) 


n  +  1 


,  2(2n  +  1)a2 

(>  -OV - 


n(n  +  1 )T 


(26) 


For  the  exterior  region,  we  have 


/  (Vw2)  •  dc  -  8*a  l  a2  ,  (27) 

n 


and  the  dissipative  factor  for  the  exterior  region  is 
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Iha  total  dissipative  factor  is  given  by  equations  (28)  and  (27),  and  if  this 
is  used  in  Lagrange's  equations. 


d  3L _ 3L_  _  _  3R 

dt  •  3a  "  • 

3a  n  3a 
n  n 


with  L  »  T  -  U#  and  R  ■  R^  +  Rq,  then  we  have 


(2n  +  1  -  nr  1 


4mr(n  -  1)(n  ♦  2)a 


n(n  +  nr 


„  ,-n  +  2^  * 

8vn  a( — — rj  a 
vn  +  1 J  n 


If  we  as sune  a  ties  dependence  of  the  form  e  in  equation  (30),  we  find 
that  the  decay  time  (imaginary  part  of  to)  is  given  by 


n  (2n  +  1)D„  ' 

Xl 


where 


>n  "  ni^n2  “  1  “  n(n  +  2)0  -  rnH  +  (2n  +  1)(1  -  T^)} 


+  n  n(n  +  2)  T 
o  n 


The  resonant  frequency,  u  (the  real  part  of  u),  is  given  by 


a> 
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2  2  (  /In  2 

“n  ■  K  ’  ”o  ■  °)  -  (r] 

n 


1  2 

for  —  <  a)  fn,  *  n  *  0)  and  w  fn.  *  n  -  0]  is  given  by  equation  (18). 
T  n'i  o  J  nv  l  o  ' 


If  in  equations  (31)  and  (32)  we  let  n£ 
given  in  Morrison  et  al,5  or 


p  *  0,  we  recover  the  result 
o 


n^n  -  1) 


>(n  -  1  )(n  +  2) 
3 


where  the  latter  result  is  that  of  Rayleigh4  and  the  decay  tine  is  the  result 

of  Lamb.6  Because  of  the  large  number  of  parameters  in  the  expressions 

for  <■)  (equation  (32))  and  the  decay  time,  t  (equation  (31)),  no  general 
n  n 

evaluation  of  the  results  is  possible.  Hence,  it  might  be  better  to  use 

equations  (31)  and  (32)  to  calculate  t  and  u  and  compare  the  results  with 

n  n 

data  taken  either  in  clouds  or  in  a  controlled  laboratory  experiment. 
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